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1.  Introduction 


The  equation  of  state  (EOS)  is  one  of  the  most  important  notions  in  thennodynamics.  A  vast 
variety  of  the  EOS  has  been  suggested  during  several  centuries  of  developing  thermodynamics. 
Various  applications  in  dealing  with  high  pressure  can  be  found  in  references  1  and  2.  During 
those  centuries,  the  very  concept  of  the  EOS  itself  was  significantly  changed,  modified,  and 
specialized  by  the  greatest  minds  of  physics,  engineering,  and  mathematics.  Because 
thennodynamics  is  a  main  interest  for  various  disciplines  and  communities,  different  meanings 
are  often  assigned  to  the  same  terms.  In  the  monograph  in  reference  3,  the  author  discusses  this 
semi-linguistic  Babel  Tower  phenomenon  regarding  the  term  “chemical  potential.”  Not 
surprisingly,  no  full  terminological  consensus  exists  between  the  authors  of  the  valuable  reviews 
(4-8).  Therefore,  in  order  to  avoid  possible  misinterpretation,  we  briefly  describe  our  own 
interpretation  of  the  existing  terminology  regarding  the  EOS. 

The  notion  of  the  EOS  appeared  well  before  the  discovery  of  the  First  Law  of  Thermodynamics. 
At  that  time,  the  EOS  was  the  function  linking  pressure  P  ,  volume  V ,  and  empirical  temperature 
r—  P  =  P(V,t)  of  a  gas  or  liquid.  Today,  many  researchers  call  this  a  thermal  EOS. 

The  discovery  of  the  First  Law  of  Thermodynamics  introduced  the  concept  of  the  internal  energy 
E  .  This  was  treated  as  a  function  of  the  volume  and  empirical  temperature  E  =  E(V,  r)  and 

often  called  the  caloric  EOS.  Functions  P(V,r )  and  E(V,t )  are  deeply  connected  with  each 
other  and  cannot  be  assigned  independently  of  each  other.  Otherwise,  the  First  Law  of 
Thermodynamics  will  be  violated.  Taken  separately,  none  of  the  functions  P(y  P)  and  E(V,t) 
are  sufficient  for  describing  all  thermodynamic  properties  of  substances. 

Later,  the  discovery  and  acceptance  of  the  Second  Law  of  Thennodynamics  entailed  introduction 
of  the  absolute  temperature  T  and  entropy  S  .  This  law  reveals  an  even  deeper  connection 
between  the  functions  of  thennal  and  caloric  EOS  than  the  First  Law  did.  Moreover,  the  Second 
Law  suggests  treating  the  internal  energy  not  as  the  function  of  the  variables  (V,r)  but  as  the 

function  of  the  pair  ( E,S) :  E  =  E(V,S) .  The  most  meticulous  theorists  call  this  the  complete 

EOS  (see  Menikoff  [3]).  Here,  “complete”  emphasizes  that  this  function  taken  alone  contains  as 
much  infonnation  as  the  two  functions  P(V,rj  and  E(V,r  )  together.  At  last,  when  the 

complete  EOS  in  the  form  E  =  E(V,S)  is  resolved  with  respect  to  the  entropy,  we  arrive  at  the 
following  EOS:  S  =  S(V,E)  .  The  last  one  is  coined  as  the  complete  EOS  for  hydrocode.  In 
this  very  sense,  the  tenn  complete  EOS  for  hydrocode  is  used  in  this  report. 

The  logical  foundations  of  classical  thermodynamics  remain  under  much  criticism — this  is  a 
natural  process  emphasizing  the  importance  of  thermodynamics.  Fortunately,  this  rarely  leads  to 
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revising  the  key  universal  identities  valid  for  complete  EOS.  Practitioners  rarely  experience 
dissatisfaction  with  the  canonical  identities  for  the  EOS.  However,  a  variety  of  novel  EOS  for 
novel  substances  is  needed.  Because  of  the  technological  and  engineering  advances,  even  for 
traditional  substances,  practitioners  demand  the  extension  of  the  earlier  accepted  EOS  for 
considerably  wider  domains  of  the  parameter  space.  Also,  there  is  an  opposite  demand  of 
developing  various  simplified  EOS  suited  for  a  specific  engineering  situation. 

From  a  physics  and  applied  mathematics  standpoint,  the  EOS  appear  to  be  simple  and 
conceptually  transparent  in  two  polar  cases.  First,  the  simplicity  is  achieved  in  the  “big  picture” 
approach  of  the  classical  thermodynamics  when  the  most  important  features  only  are  explicitly 
taken  into  account,  whereas  the  variety  of  minor  features  deliberately  remain  unspecified. 
Second,  the  simplicity  appears  in  the  opposite  case  of  extremely  specialized  models  (e.g.,  the 
perfect  gas  model).  Besides  being  useful  in  understanding  physical  phenomena,  these  cases  also 
allow  deep  theoretical  analysis. 

It  is  now  universally  accepted  that  no  single  EOS  can  work  for  all  purposes,  all  substances,  and 
all  ranges  of  thermodynamic  parameters.  Even  if  such  a  universal  model  were  developed,  it 
would  only  be  useful  for  blind  computer  implementations  and  not  for  further  theoretical  analysis 
or  intuitive  comprehension.  Theoretical  analysis,  however,  is  very  important,  even  in  the 
computer  era. 

We  suggest  a  novel  model  that  is  rather  simple  and  conceptually  transparent.  It  possesses  some 
features  of  the  big  and  small  pictures  previously  mentioned.  It  has  infinitely  many  parametric 
degrees  of  freedom.  However,  the  infinitely  many  degrees  of  freedom  appear  not  in  the  form  of 
infinitely  many  unrelated  parameters  but  in  the  form  of  two  functional  degrees  of  freedom.  This 
model  is  complete  in  the  thermodynamic  sense  and  includes  some  classical  models  as  special 
cases.  At  the  same  time,  it  appears  to  be  convenient  for  converting  the  shock-wave  locus  or  the 
static  adiabata  locus  into  the  EOS.  Explicit  formulas  of  this  conversion  are  given  next. 

We  call  our  models  operational  EOS’s.  The  term  operational  is  used  because  of  two  reasons. 
First,  the  suggested  EOS’s  mostly  have  the  fonn  of  various  operators  with  respect  to  the 
experimental  data  used  in  their  definitions.  Second,  they  have  direct  relation  to  the  Bridgman 
methodology  of  operationalism.  We  will  discuss  these  methodological  issues  in  forthcoming 
reports  after  demonstrating  several  important  examples. 


2.  The  Suggested  Model 


The  suggested  model  is  not  as  general  as  the  most  general  models  of  classical  thermodynamics. 
In  fact,  we  suggest  a  class  of  the  EOS  for  the  hydrocode  based  on  the  following  two  simple 
constraints: 
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1 .  The  heat  capacity  at  constant  volume  Cv  is  a  function  of  the  internal  energy  density  E 
only: 

Cr  =  x(E).  (1) 

2.  The  nondimensional  combination  PV  /  RT  is  a  function  of  the  volume  V  only: 


PV 

RT 


Ab)> 


(2) 


where  R  is  the  universal  thermodynamic  constant. 

Some  classical  EOS’s  are  automatically  included  in  this  category.  For  instance,  the  model  of  the 
monoatomic  perfect  gas  belongs  to  this  category  with  the  following  choice  of  the  parametric 
functions: 


*(£)  =  y.A(r)  =  l. 


(3) 


Another  example  is  the  monoatomic  van  der  Waals  gas  with  the  vanishing  correction  coefficient 
for  the  internal  gas  pressure.  This  classical  model  corresponds  to  the  following  choice  of  the 
functions: 


x(e)  =  \r,  Mr) 


V 

V-b 


(4) 


In  view  of  what  was  discussed  in  the  Introduction,  even  the  most  innocent-looking  assumptions 
can  be  thermodynamically  inconsistent.  In  other  words,  there  may  not  exist  any  complete  EOS 
for  which  those  assumptions  can  be  deduced  in  a  thermodynamically  consistent  way.  It  can  be 
demonstrated,  however,  that  the  assumptions  1  and  2  are  thermodynamically  consistent  for  any 
of  the  two  meaningful  functions  %(E)  and  A(F) .  Moreover,  such  a  potential  can  be  presented 

in  the  following  explicit  form: 


S(V,E)-S‘  = 


RA(m) 


</£ 


m 


Ih{Z) 


(5) 


where  A  is  a  constant.  The  function  H  has  the  physical  meaning  of  the  absolute  temperature. 

A  simple  verification  of  thermodynamic  consistency  and  validity  of  equations  1,  2,  and  5  is 
presented  in  section  2.2. 

For  the  gas  model  described  by  equation  4,  the  EOS  implies  the  following  complete  EOS  for 
hydrocode: 
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(6) 


S  (  F,  is )  =  i?  In 


V-b 
V"  -b 


3  R  E-E.  +  CA 

+  ~CV  In - * - — 

2  CrA 


2.1  Recovery  of  the  Functions  %(E)  and  A(F)  From  Measurements 

The  class  of  substances  described  by  the  A/  -  EOS  is  much  wider  than  the  classical  models.  In 
fact,  this  class  contains  two  arbitrary  functions  which  can  be  specified  based  on  the  available  sets 
of  experimental  data. 


2.1.1  The  Hugoniot  Adiabata  Data 

Given  the  Hugoniot  adiabata  locus  EH  ( V) ,  we  get  the  following  Rankine -Hugoniot  relationship 
(8-10): 


2EH(V)  +  P(V,EH(V))(V-K)-2E.+P.(V-K)=0,  (7) 


which  implies  the  following  formula  of  A(F) : 

Air)  2Eh(V)-2E.-P.(K-V)  V 
1  '  RH(Eh(V ))  K-V 


(8) 


and  the  following  complete  EOS  for  hydrocode: 

1  p  .  2EH(m)-2E,-P.(V.-m) 


S(V,E)  =  S*+  [  ds — — —  +  f  dm 
K  ’  i*  H(e)  J 


(V.-m)H{EH(m)) 


(9) 


In  the  case  of  vanishing  ambient  pressure  P,  =  0  and  when  the  heat  capacity  Cv  is  constant,  we 
arrive  at  the  following  simple  EOS: 


S(V,E)  =  S,  +  Cv  \nE  E*  +  CvT 


CVT 


J  dur 


2EH(tu)-2Et 


(y.-m)\EH{G>)-E.  +  CVT* ]  ’ 


(10) 


where  T  is  the  ambient  temperature. 

Equations  9  and  10  show  that  inversion  of  the  experimental  measurements  of  the  heat  capacity 
Cv  (£■)  and  the  Hugoniot  locus  EH  (V )  into  the  complete  EOS  for  hydrocode  S(V,E)  is  fulfilled 

with  the  help  of  the  integral  operator.  The  operator  has  singularity  in  the  vicinity  of  the  ambient 
configuration. 


2.1.2  The  Quasistatic  Adiabata  Data 

The  kernel  of  the  integral  operator  in  the  recovery  formula  has  a  singularity  near  the  ambient 
state,  i.e.,  for  the  weak  shock  waves.  Therefore,  in  the  vicinity  of  the  ambient  state,  some  other 
approaches  to  recovering  the  EOS  should  be  used. 
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The  quasistatic  adiabata  is  given  by  equation  E  =  E*S(V).  Using  equation  5,  we  get  the 
following  inversion: 

E  1 

S(V,E)=  [  d% - ? - +  5*. 

E's(V)  T*  +  f  ^ 

l.x{v) 


(11) 


In  the  case  of  constant  heat  capacity  at  fixed  volume,  equation  1 1  reads  as  follows: 


S(V,E)  -  S* 


C  In  E~E*  +  Cvr 
v  El(V)-E«  +  Cvr ‘ 


(12) 


Equations  1 1  and  12  show  that  inversion  of  the  experimental  measurements  of  the  heat  capacity 
CV(E )  and  the  static  loading  locus  Es  (V)  into  the  complete  EOS  for  hydrocode  S(V,E  )  is 

fulfilled  with  the  help  of  other  operators,  as  compared  with  equations  9  or  10.  The  operators 
remain  integral  in  the  case  of  nonconstant  heat  capacity.  For  the  case  of  constant  heat  capacity, 
the  operator  becomes  functional  but  transcendental.  Singularities  in  this  case  appear  when  the 
ambient  absolute  temperature  is  equal  to  O  K  as  many  reserchers  prefer.  Of  course,  the 
assumption  of  constant  heat  capacity  does  not  work  in  the  vicinity  of  O  K . 

For  explicit  formulas  for  other  types  of  experimental  data,  see  Grinfeld  (11). 


2.2  Verification  of  Equations  1  and  2  for  the  EOS  (5) 

Combining  equation  3  with  the  thermodynamic  identities 

dS(E,V)_  1  dS(E,V)_P 
dE  ~  T  ’  dV  ~  T  ’ 


we  get 


(13) 


(14) 


Thus,  the  absolute  temperature  for  the  model  under  discussion  depends  upon  the  internal  energy 
only,  and  constant  A  has  the  physical  meaning  of  the  absolute  temperature  at  E  =  E * . 


Equation  14  proves  equation  2.  Combining  equation  14  with  the  thermodynamic  identity 


cv  = 


dT(V,E ) 


dE 


(15) 


we  arrive  at  equation  2: 


CV=Z(E)  =  CV(E). 


(16) 
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3.  Discussion  and  Conclusion 


We  suggest  a  novel  class  of  simple  EOS  for  hydrocode,  depending  on  two  functions  of  one 
variable  each.  The  suggested  A/  -  models  generalize  some  classical  models  but  also  allow  a 

wider  set  of  degrees  of  freedom.  Functions  A(V)  and  /(T)  should  be  recovered  by  comparing 

available  experimental  data. 

This  is  not  the  most  general  form  of  the  complete  EOS.  Some  data,  when  available,  cannot  be 
considered  in  the  suggested  framework.  A  much  more  general  class  of  the  EOS  is  suggested  in  a 
future  report  (11).  On  the  other  hand,  the  suggested  model  has  its  own  advantages.  First,  it 
allows  further  analytical  consideration.  Second,  it  can  be  recommended  in  cases  where  the  scope 
of  data  is  relatively  narrow. 

We  established  the  explicit  formulas  for  the  recovery  of  the  complete  EOS  for  hydrocode  from 
experimental  data  for  the  heat  capacity  at  constant  volume  Cv  and  the  loci  of  the  Hugoniot  or 

static  adiabata.  The  recovery  formulas  have  the  form  of  integral  operators  with  certain 
singularities.  Those  formulas  can  be  relatively  easily  implemented  into  computer  code. 

Vigilance  is  necessary  in  handling  the  singularities. 

The  suggested  EOS  is  complete  in  the  sense  that  it  automatically  satisfies  the  First  and  Second 
Laws  of  Thermodynamics. 
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LOS  ALAMOS  NATIONAL  LAB 
SR  GREENFIELD  MS  J565 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
R  MENIKOFF  MS  B214 
LOS  ALAMOS  NM  87545 


LAWRENCE  LIVERMORE  NATL  LABS 
L  BENEDICT 
M  A  BARRIOS 
A  TEWELDEBERHAN 
7000  EAST  AVE 
LIVERMORE  CA  94550 


LOS  ALAMOS  NATIONAL  LAB 
B  J  JENSEN  MS  B952 
LOS  ALAMOS  NM  87545 


NAVAL  SURFACE  WARFARE  CTR 
RESEARCH  AND  TECHLGY  DEPT 
G  T  SUTHERLAND 
INDIAN  HEAD  MD  20640 


LOS  ALAMOS  NATIONAL  LAB 

J  N  JOHNSON 

LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
B  PLOHR  MS  B213 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
L  SMILOWITZ  MS  P915 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
D  PRESTON 

LOS  ALAMOS  NM  87545 

S ANDIA  NATIONAL  LABS 
J  H  CARPENTER  MS  1322 
ALBUQUERQUE  NM  87185 

S ANDIA  NATIONAL  LABS 
M  DESJARLAIS  MS  1189 
TR  MATTSON  MS  1189 
ALBUQUERQUE  NM  87185 

S ANDIA  NATIONAL  LABS 
RJ  MAGYAR  MS  1322 
ALBUQUERQUE  NM  87185 


1  APPLIED  RESEARCH  ASSOC 
D GRADY 

4300  SAN  MATEO  BLVD  A-220 
ALBUQUERQUE  NM  871 10 

1  US  AIR  FORCE  RESEARCH  LAB 
MUNITIONS  DIRCTRT 

Y  HORIE 

EGLIN  AFB  FL  32542 

2  THE  JOHNS  HOPKINS  UNIV 
DEPT  OF  MECH  ENGRG 

K  T  RAMESH 
T  W  WRIGHT 
232  LATROBE  HALL 
3400  N  CHARLES  ST 
BALTIMORE  MD  21218 

2  UNIVERSITY  AT  STONY  BROOK 
J  GLIMM 
R  SAMULYAK 

STONY  BROOK  NY  11794-3600 

2  UNIVERSITY  OF  SOUTH  FLORIDA 
DEPT  OF  PHYSICS 
I  OLEINIK 

V  ZHAKHOVSKY 
TAMPA  FL  33620 
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1  UNIV  OF  UTAH 

DEPT  OF  MECH  ENGRG 
R  BRANNON 

50  S  CENTRAL  CAMPUS  DR 
SALT  LAKE  CITY  UT  841 12 

1  UNIV  OF  CALIFORNIA 

DEPT  OF  MECH  &  ARSPC  ENGRG 
V  NESTERENKO 
9500  GILMAN  DR 
LA  JOLLA  CA  92093-0411 

1  UNIV  OF  MD 

DEPT  OF  MECHANICAL  ENGRG 
R  ARMSTRONG 
COLLEGE  PARK  MD  20742 

1  UNIV  OF  COLORADO  AT  BOULDER 
DEPT  OF  CIVIL  ENGRG  & 
ARCHITECTURAL  ENGRG 
R  A  REGUEIRO 
428  UCB  ECOT  441 
1 1 1 1  ENGRG  DR 
BOULDER  CO  80309-0428 

1  ROBERT  R  MCCORMICK  SCHOOL 
OF  ENGRG  &  APPLIED  SCI 
H  ESPINOSA 
633  CLARK  ST 
EVANSTON  IL  60208 

1  UNIV  OF  MISSOURI 
DEPT  OF  CHEMISTRY 
T  SEWELL 

COLUMBIA  MO  652 1 1  -7600 
1  UNIV  OF  MD 

DEPT  OF  MATERIALS  SCI  &  ENGRG 
M  KUKLJA 

COLLEGE  PARK  MD  20742 

1  STANFORD  UNIV 

DEPT  OF  AERONAUTICS  & 
ASTRONAUTICS 
R  CHRISTENSEN 
DURAND  BLDG  496 
LOMITA  MALL 
STANFORD  CA  947305-4035 

1  US  ARMY  RESEARCH  OFFICE 
RDRL  RO 
D  SCUTRUD 
BLDG  4300 
DURHAM  NC  27703 


1  US  ARMY  RESEARCH  OFFICE 
RDRL  ROI 
B  WEST 
BLDG  4300 
DURHAM  NC  27703 

1  US  ARMY  RESEARCH  OFFICE 
RDRL  ROI  M 

J  MYERS 
BLDG  4300 
DURHAM  NC  27703 

2  US  ARMY  RESEARCH  OFFICE 
RDRL  ROE  M 

J  PRATER 
D  STEPP 
BLDG  4300 
DURHAM  NC  27703 

1  US  ARMY  RESEARCH  OFFICE 
RDRL  ROE  N 
R  ANTHINIEN 
BLDG  4300 
DURHAM  NC  27703 

3  US  ARMY  RESEARCH  OFFICE 
J  LAVERY 

D  MANN 
S  MATHAUDHU 
BLDG  4300 
DURHAM  NC  27703 

ABERDEEN  PROVING  GROUND 

72  DIR  USARL 
RDRL  WM 
B  FORCH 
J  MCCAULEY 
SKARNA 
P  PLOSTINS 
RDRL  WMP  B 
R  BECKER 
S  BILYK 
D  CASEM 
J  CLAYTON 
D  DANDEKAR 
J  FITZPATRICK 
M  GRINFELD  (20  CPS) 

C  HOPPEL 
R  KRAFT 
B  LEAVY 
D  POWELL 
M  RAFTENBERG 
M  SCHEIDLER 
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T  WEERASOORIY A 
C  WILLIAMS 
RDRL  WMP  C 
T  BJERKE 
S  SEGLETES 
W  WALTERS 
RDRL  WML  B 
I BATYREV 
J  BRENNAN 
S  IZVEKOV 
B  RICE 
RDRL  CIH  C 
P CHUNG 
J  CAZAMIAS 
J  KNAP 
RDRL  WMM  G 
I  ANDZELM 
F  BEYER 
J  LENHART 
B  RIDERSPACHER 
RDRL  WMP  E 
W  GOOCH 
B  LOVE 

C  KRAUTHAUSER 
RDRL  WMM  E 
J  ADAMS 
L  KECSKES 
J  LASALVIA 
P  PATEL 
J  SWAB 
M  WILL-COLE 
RDRL  WMM  B 
B  CHEESEMAN 
C  FOUNTZOULAS 
G  GAZONAS 
RDRL  SER  L 
W  NOTHWANG 
RDRL  WML  H 
M  FERMEN-COKER 
RDRL  WMP 
B  BURNS 
S  SCHOENFELD 
RDRL  WMP  D 
H  MEYER 
RDRL  WMP  G 
R  BANTON 
S  KUKUCK 
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1  CAVENDISH  LAB 

MOTT  BLDG  RM316B 
JJ  THOMSON  AVE 
M  CHAUDHRI 
CAMBRIDGE  CB3  OHE 
UNITED  KINGDOM 

1  LOUGHBOROUGH  UNIV 

WOLFSON  SCHOOL  OF  MECH  & 
MANUFACTURING  ENGRG 
LOUGHBOROUGH 
V  SILBERSCHMIDT 
LE11  3TU 

UNITED  KINGDOM 

1  UNIVERSITY  OF  LIVERPOOL 
DEPT  OF  ENGRG 
BROWNLOW  HILL 
LIVERPOOL 
H  OUYANG 
L69  3GH 

UNITED  KINGDOM 

4  AWE  ALDERMARSTON  READING 
BERKSHIRE  RG7  4PR 
N  BOURNE 
J  C  F  MILLETT 
G  A  COX 
C  M  ROBINSON 
UNITED  KINGDOM 

1  DEFENCE  SCI  &  TECHLGY 
ORGANISATION 
WEAPONS  SYSTEMS  DIV 
A  RESNYANSKY 
EDINBURGH  SA5111 
AUSTRALIA 

1  CSIRO  EXPLORATION  &  MINING 

PO  BOX  883 
KENMORE  QLD  4069 
F  D  STACEY 
AUSTRALIA 

1  MOSCOW  STATE  UNIV 

INSTITUTE  OF  MECHANICS 
MICHURIN SKI(E  PR  1 
S  S  GRIGORYAN 
MOSCOW  117192 
RUSSIA 


1  STEKLOV  MATHEMATICAL 
INSTITUTE 
GUBKINA  AV  8 
A  G  KULIKOV SKII 
MOSCOW  117966 
RUSSIA 

1  LANDAU  INSTITUTE  FOR 
THEORETICAL  PHYSICS 
N  INOGAMOV 

RAS  CHERNOGOLOVKA  142432 
RUSSIA 

5  JOINT  INSTITUTE  FOR  HIGH 
TEMPERATURE  PHYSICS 
RUSSIAN  ACADEMY  OF  SCIENCES 

V  FORTOV 

V  GRYAZNOV 
GKANEL 

M  LOMONOSOV 

S  RAZORENOV 

JIHT  RAS  MOSCOW  125412 

RUSSIA 

1  INSTITUTE  OF  CONTINUOUS  MEDIA 
MECHANICS 
URAL  BRANCH  OF  RAS 
O  B  NAIM  ARK 
1  ACAD  KOROLEV  STR 
PERM  614013 
RUSSIA 

1  BEN  GURION  UNIV 

DEPT  OF  MECHANICAL  ENGRG 

E  B  ZARETSKY 

PO  BOX  653 

BEER  SHEVA  84105 

ISRAEL 

1  RAFAEL 
PO  BOX  2250 

V  PARTOM 
HAIFA  3 1021 
ISRAEL 
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